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Q.J Abstract: We analyze the validity of BCFW recursion relations for currents of n — 2 gluons and two 
massive quarks, where one of the quarks is off shell and the remaining particles are on shell. These currents 
are gauge-dependent and can be used as ingredients in the unitarity-based approach to computing one-loop 
amplitudes. The validity of BCFW recursion relations is well known to depend on the so-called boundary 
l/") . behavior of the currents as the momentum shift parameter goes to infinity. With off-shell currents, a 

m 



new potential problem arises, namely unphysical poles that depend on the choice of gauge. We identify 
, conditions under which boundary terms are absent and unphysical poles are avoided, so that there is a 
natural recursion relation. In particular, we are able to choose a gauge in which we construct a valid shift 
for currents with at least n — 3 gluons of the same helicity. We derive an analytic formula in the case where 
all gluons have the same helicity. As by-products, we prove the vanishing boundary behavior of general 
• rH off-shell objects in Feynman gauge, and we find a compact generalization of Berends-Giele gluon currents 
^ . with a generic reference spinor. 
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1. Introduction 



Scattering amplitudes in gauge theories rapidly become difficult to compute as the number of external 
legs increases. The difficulty is encountered especially when seeking analytic expressions. However, recent 
advances involving on-shell techniques have enabled the discovery of many new formulas for amplitudes. 
Notably, the BCFW construction generates tree-level amplitudes efficiently and compactly through recur- 
sion relations, via the Cauchy residue theorem. 

The phrase "on-shell technique" refers to methods of computing scattering amplitudes in which certain 
propagators are taken to their on-shell limits, and amplitudes are then constructed from knowledge of the 
corresponding factorization properties. The key is that these factorization limits are calculated in terms 
of amplitudes of lower complexity, i.e. fewer legs or loops. 

It is interesting to study amplitudes (or rather, currents) in which one or more legs is continued off shell, 
since they carry even more information than on-shell amplitudes. For example, the Berends-Giele recursion 
relation among gluon currents in Yang-Mills theory [1] is not only computationally powerful for numerical 
results, but was also the crucial stepping stone to establishing the first formulas for gluon amplitudes with 
arbitrary numbers of legs, in certain helicity configurations [1,2]. It is still possible to consider the limits 
in which internal propagators go on-shell and apply the BCFW construction to find recursion relations [3] . 
Compared to the recursion relations for on-shell amplitudes, the ones for currents require committing to a 
gauge choice, and summing over all internal polarization states, including unphysical polarizations. 

In this paper, we seek compact analytic forms for currents of n — 2 gluons and two massive quarks, 
where one of the quarks is off shell and the remaining particles are on shell. These currents are key 
ingredients of an on-shell method of computing 1-loop amplitudes with external massive fermions [4] . Such 
amplitudes are of particular interest in the context of LHC searches for new physics, where production of 
top quarks plays a large role in both signals and backgrounds. Massive fermion currents can be computed 
with the off-shell Berends-Giele recursion [1]. In [5], this was used to give a compact result in the case 
where all gluons have the same helicity, with a particular gauge choice relative to the massive spinors. 

We study the validity of the BCFW construction [6, 7] for these massive fermion currents. The con- 
struction begins by shifting the momenta of a pair of on-shell external legs by +zq and — zq respectively, 
where z is a complex variable and q is obtained by requiring that both legs remain on shell after the shift. 
Then, the residue theorem produces a recursion relation from poles in z taking values where propagators 
go on shell. The construction breaks down if there are poles from other sources. In Yang-Mills theory, the 
only other possible source is a "boundary term," from a nonvanishing limit as z is taken to infinity. 1 For 
off-shell currents, there is another problematic source of poles, which we call "unphysical poles." They are 
due to the dependence on gauge choice, and they spoil the recursion relation, since we have no information 
about how to calculate their residues independently. 

x If the theory is sufficiently well understood, it is possible to include a boundary term explicitly at each step of the 
recursion [3,8-11]. 
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We identify conditions under which the boundary terms and unphysical poles vanish for massive 
fermion currents, so that the BCFW construction produces a recursion relation. We then proceed to solve 
the recursion relation in the particular case where all gluons have equal helicities. Compared to the more 
compact result of [5], our formula also requires all gluons to use the same reference spinor but preserves 
the genericity of its value. 

Our analysis of boundary terms is based on grouping Feynman diagrams conveniently and applying 
the Ward identity and inductive arguments. The argument establishes the absence of boundary terms 
for general off-shell objects in Feynman gauge, provided that there are two on-shell gluons available to 
construct the momentum shift. 

In our study of unphysical poles, we use off-shell gluon currents of the type originally derived by 
Berends and Giele [1]. We are motivated to generalize the currents in which one gluon has opposite 
helicity to all the others, by taking its reference spinor to be arbitrary. When the opposite-helicity gluon is 
color- adjacent to the off-shell leg, we find a very compact form for the current. When it is centrally located 
among the other gluons, we prove that the current is, in fact, independent of the arbitrary reference spinor. 

This paper is organized as follows. The remainder of this introduction contains details of our con- 
ventions and notation. Section 2 reviews the BCFW construction of recursion relations in the context of 
massive fermion currents and describes the origin of boundary terms and unphysical poles. Section 3 derives 
sufficient conditions for good boundary behavior; Section 4 derives sufficient conditions for the absence of 
unphysical poles. Section 5 presents sample results from our recursion relations. In particular, we find a 
closed form for n-point currents in which all gluons have the same helicity. Section 6 is a summary with 
proposals for future work. Appendix A works out a technical point in the proof of Section 3, namely the 
use of Ward identities in Feynman gauge. Appendix B presents a generalization of Berends-Giele currents 
with one gluon of opposite helicity, in which the choice of reference spinor is relaxed. Appendix C contains 
a fully non-recursive formula as an alternative to the n-point current given in Section 5, and outlines its 
derivation. 

For reference, analytic formulas for on-shell amplitudes of gluons with massive quarks may be found 
in [12-20], all of which use BCFW recursion relations, sometimes in combination with SUSY Ward identities 
or Berends-Giele recursion. 

1.1 Conventions and notation 

Momenta of gluons are directed outward, while momenta of fermions are directed inward. We will be 
considering color-ordered amplitudes and off-shell currents with one massive fermion line, for example, 
iJ (1q, 2q, S g , A g , . . . ,rig^, where the star means that the indicated leg is considered off-shell, while the 
remaining legs are on-shell. We do not include the propagator for the off-shell leg in our definition. For this 
current, the quark line has its arrow pointing from leg 2 to leg 1. When the quark line matrices are read 
against the arrow, then the gluon indices are contracted in reverse numerical order. In slightly different 



- 3 - 



notation, we can write 

iJ fe, 2 Q , 3 h *,4 h g \..., n h A = \n h - . . . 4^ |2), (1.1) 



Q' <V' "9 ' 9 ' • • • ' "g 

where the round bracket |2) can be equal to either |2) or |2], depending on its spin. This notation emphasizes 
the fact that the current is a spinorial object. For example, to obtain the corresponding amplitude, one 
should first put p\ on shell and then contract the current with either [1| or (1|. 
Our color-ordered Feynman rules use the following gluon vertices: 

q,H r,u 

= [SV(P - q)u + 9^{q - r)\ + g v \(r - p) M ] , (1.2) 

p,\ 
P- v 

^yC^ = 7; [2g\v9w ~ 9x P 9^u - 9x^9u P ] • (1.3) 
A P 

The polarization vector for a gluon of momentum p is, depending on helicity [21-25], 

n 1 [n p \Y\p) n 1 {n p \^\p\ 

p— /k r„ „i ' p+ /k /„ „\ ' l 1 - v 



where n p is an arbitrary but fixed "reference" momentum satisfying r? v = and either (n p p) / or 
[n p p] / 0, so that the denominator is nonzero. The null reference momenta are chosen independently for 
each gluon. The set of reference momenta is what we refer to as the gauge choice for a particular calculation, 
within the Feynman gauge used throughout the spinor-helicity formalism. Any current we construct with 
a specific gauge choice is expected to fit into a larger calculation, such as the one-loop computations of [4], 
in which all external legs are on-shell, so that ultimately, after being combined with other ingredients 
computed in the same gauge, no trace of the gauge choice remains. Therefore the reference spinors can be 
chosen to maximize computational convenience. We delay the choice as far as possible, so that convenience 
can be evaluated later in the full context of a larger calculation. 

The spinors for the massive fermions satisfy the Dirac equation. We do not require any further details 
of their definitions, so any of various conventions (e.g. [23], [26]) can be used. The massless limit is smooth. 

We use the Lorentz vector Pij to denote the sum of color-adjacent momenta in increasing cyclic order, 
between and including legs i and j. 

2. The recursive construction 

We choose to apply the momentum shift to a pair of gluons, since they are always on-shell in this context. 
The shift denoted by [kl) shifts the momenta of gluons labeled by k and I as follows: 

\k) = \k), \k] = \k]-z\i], \i) = \i) + z\k), \i] = \i], (2.i) 
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or 

Pk=Pk~ zq, Pi =Pi + zq, (2.2) 

with the complex-valued shift vector 

t = \ <fc|7l*] • (2-3) 
The BCFW construction is to apply the residue theorem on iJ{z)/z to recover the current as 

iJ(z = 0) = - ^ Res 

poles 

To obtain recursion relations among currents, we must have no pole at infinity. This is assured if J(z) 
goes to zero in the limit z — > oo, which we call good boundary behavior. 

The remaining poles can have two possible origins, due to their construction from Feynman rules: 
(1) the vanishing of a propagator, or (2) the vanishing of the denominator of a polarization vector, when 
written as in (|1.4|). 

The first type of pole, from a vanishing propagator, is familiar from the recursion relation for on-shell 
amplitudes. The corresponding residues are easy to evaluate as the product of two currents or amplitudes 
with fewer legs, since the vanishing of the propagator is an on-shell condition. 

The second type of pole will be called an unphysical pole. In an on-shell amplitude, the reference 
spinor in such a pole could be freely chosen to eliminate the z dependence in the denominator, but now 
with a quark line off-shell, we must fix all reference spinors from the start, and they play an explicit role. 
These unphysical poles are problematic, since their location has no natural physical meaning, and we have 
no independent way of computing their residues. Thus, we will find conditions on the currents and shift 
that prevent the appearance of unphysical poles. 

With good boundary behavior and no unphysical poles, there will be a recursion relation that takes 
the schematic form 

ij n (z = 0)= ^2 ^2iJ k+1 (z) iJ n -k+i(z), (2-5) 

where the iJ are currents and amplitudes with fewer legs; the hat denotes their evaluation at shifted 
momentum values; P is the momentum flowing between them which goes on-shell at the pole z; and the 
second sum is over internal helicities. Of course, the propagator acquires a mass if it is fermionic. If all of 
the off-shell legs belong to just one of the two currents on the right-hand side, then the other is replaced 
by a shifted on-shell amplitude, iM(z). 

3. Boundary behavior 

In this section, we study the behavior as z — > oo of the fermionic current iJ (^q, 2q, 3 9 , 4 g , . . . , n g ^ under 
the \kl) shift ( |2.1[) , where k and / represent any of the gluons. We will conclude that the boundary term 



- 5 - 



vanishes in the helicity cases (k~,l£), (k~,l~) and thus as well, for a generic gauge choice. In fact, 

our argument is much more general, since none of the unshifted gluons need to be on shell. Moreover, the 
number of fermion lines is not crucial either. 



3.1 Choice of shift: helicities and polarizations 

Consider the superficial boundary behavior of individual Feynman diagrams, following the flow of the 
additional momentum zq. Without the polarization vectors, the diagrams where the zq momentum goes 
only through 3-gluon vertices and gluon propagators behaves the worst — as O(z). If zq runs through a 
4-gluon vertex or through a fermion line, then such a diagram already behaves as 0(1) or better. Now, if 
we contract the vector indices for the [kl) -shifted gluons with their polarization vectors, 

^ 1 [n k mk) p 1 (n k \^\k] - z(n k \^\l] 

k ~ V2 [n k k] - z[n k l] k+ y/2 (n k k) 

~u l [ ni \Y\l) + z[m\Y\k) l (ni\r\q 

y/2 [nil] l+ V2(nil)+z(mk)' 

we see that for a generic gauge choice the off-shell current superficially has O (~) behavior in the (k~ , Z+) 
case, O(z) in the (k~,l~) and cases, and 0(z 3 ) in the (k^,l~) case. Note that this behavior can 

be altered by special gauge choices, i.e. if = I or ni = k, then and lose their z dependence 
instead of being O (~). 



The helicity case (k g ,/+) is thus safe automatically, for a generic gauge choice. In Subsection 3A, we 
will discuss boundary behavior for a special gauge choice that will be needed in the following sections. 

In the remainder of this section, we will prove that in a generic gauge (where n k ^ I), the off-shell 
current with helicities (k~,l~) also vanishes at infinity, at least as O Q). To do that, let us multiply this 
current by the z-independent factor —y/2[nil], so that only the numerator, [nz^lZ), remains contracted 
with the Z-th gluon's Lorentz index v. The resulting expression depends only linearly on [n/|, which is a 
2-dimensional massless spinor and thus can be expressed as a linear combination of any two independent 
spinors of the same kind: 

[m\ = a[l\+l3[n k \. (3.2) 

Therefore it is enough to show that we get O Q) behavior for the two special cases n\ = I and n/ = n k . 
Let us examine them one by one. 

3.2 Like-helicity shift, first term 

The first term in ( |3.2[) yields [Z|7"|Z) = 2Z V , making it possible to use the Ward identity, which diagram- 
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matically can be expressed as follows: 




(3.3) 



n + I 



Each diagram on the right-hand side of (|3.3|) is supposed to have an appropriate gauge group generator 
contracted with the leg to which the momentum I is added, but that is irrelevant for the discussion of z — > oo 
behavior. For any leg that is initially on shell, such as the k-ih gluon, the corresponding right-hand-side 
term should naturally vanish, because the resulting leg would go off shell and thus would be left out when 
extracting the on-shell pole residue according to the standard LSZ procedure. 

Now consider what the Feynman rules tell us about the diagrams on the right-hand side. In case there 
are off-shell gluons, if a diagram has the zq momentum going from the fc-th gluon to another through 

3- gluon vertices, it must now behave no worse than O (-)! Indeed, it still has ~ O Q) on the k-th leg, 
and in addition to that the gluon propagator on the off-shell leg is now O Q) as well. If zq runs through a 

4- gluon vertex or through a fermion line and ends up still on a gluon leg, then such a diagram will behave 
at most as O (^)- A new ingredient here is the diagrams that have zq momentum flowing through 3-point 
vertices to an off-shell fermion leg, but it is easy to see that they will also behave like O (-) or better. In 
sum, applying the Ward identity in this way reduces the maximal superficial power of z at infinity by two. 

There is one technical caveat about this argument: strictly speaking, the Ward identity ( |3.3| ) is valid 
for a ghostless gauge, whereas in Feynman gauge it is necessary to introduce some extra terms on the 
right-hand side. We address this issue carefully in Appendix [A] and find that the argument still holds in 
Feynman gauge. 



3.3 Like-helicity shift, second term 

Now we consider the case ni = = n. It turns out to be possible to deduce some interesting facts about 
the boundary behavior of an off-shell current simply from the Feynman rules. 

3.3.1 Gluon trees, leading power of z 



We start with the leading 0{z) diagrams, in which the zq momentum flows only through 3-gluon vertices 



-7- 



and gluon propagators, which thus behave as O(z) at most. Let us look closely at the part of such a 
diagram that is directly adjacent to the zq momentum flow, i.e. just a gluon tree with all but the k-th and 
Z-th legs off-shell and their propagators amputated, keeping in mind that any of the off-shell legs can be 
extended by any sort of ^-independent tree, including a fermion line. A tree with n legs will have (n — 2) 



A2 A 3 




A 5 

Figure 1: A generic gluon tree diagram with only 3-point vertices. 

vertices ~ O(z), (n — 3) internal propagators ~ O and (n — 2) free indices. The highest power of z 
will be accumulated if we pick up zq x from each vertex, a zfn^j/c) term from the l-th gluon's polarization 
vector, and another [n|7 M |/c) coming from the k-th gluon: n vectors in total. Apart from that, vertices and 
propagator numerators can only offer various combinations of metric tensors, and the fact that there are 
only (n — 2) free indices means that at least one contraction will take place among those vectors. But any 
such contraction eliminates a power of z, since 

[n\j v \k) gix =0, 
q 2 = 0. 

So the leading O(z) term vanishes and we are left only with 0(1) at most. 
3.3.2 Fermion line insertion, leading power of z 

Similarly, the leading 0(1) term vanishes for the diagrams in which zq flows through the fermion line. To 
see this, consider once again only the terms directly adjacent to the zq momentum flow, i.e. the relevant 
part of the fermion line and mostly off-shell gluon trees on both of its sides. 



k A 2 A 3 i 




Figure 2: A generic diagram with zq flowing through the fermion line and only 3-point gluon vertices. 

The leading power of z in a generic diagram with n gluons will now be attained by accumulating (n — 1) 
powers of zq x from 3-gluon vertices and the numerators of fermion propagators. As before, [n|7^|A:) and 
z[n|7"|A;) will come from the k-th and l-th. gluons respectively, so in total we have (n + 1) vectors with only 
(n — 2) free indices to attribute to them. Note that the relevant part of the fermion line consists of an odd 
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number of 7-matrices and thus can always be expressed as a linear combination of eight basic matrices 
{7 M }^ =0 U {7 M 7 5 }^=o just by using the standard formula: 

tVV = 2 Aa V " 9 Xv Y + <fV + ie x ^ lpl \ (3.4) 

The free index of the 7^ or 7^7 5 can either be left free as an off-shell gluon index (leaving us with only (n— 3) 
free indices left for (n + 1) vectors) or be contracted with one of the {n + 1) vectors. So the number of free 
indices is smaller than the number of vectors at least by two. The difference with the previous case is that 
now we have not only metric tensors to do the index-contraction work, but also the totally antisymmetric 
tensor coming from (|3.4| ). So, lowering the number of free vector indices by two can be achieved by either 
dotting one vector to another, in which case we get zero just as in the gluon-tree case; or by contracting 
three vectors to one antisymmetric tensor constructing terms like ex^up • Z Q X • [^17^1^) • z [ n \l v \k), all of 
which vanish due to the fact that we have copies of only two vectors in the leading O(z) term. Thus all 
diagrams with zq momentum flow through a fermion line necessarily vanish at z — > 00. 

3.3.3 Gluon trees, next-to-leading power of z 

What is left to consider is the possible 0(1) contribution from gluon trees. To begin with, calculate the 
simplest gluon tree, i.e. a single 3-gluon vertex contracted with the two polarization vectors as given in 



(» k I 

T 

A = l — (*<fc|7AM • 2{k\l\n]) + <%A|n](J|7>] • z{k\ lv \l]) + O (3-5) 



C- (a^(k - l) x + g vX {2l + k)„ - g^(2k + Z)„) e\__ 



2z[nl]< 



''" Ai '^ (2{kl)[ln\ + 2{lk)[ln\) + O ( -) = O ( - 



2[nl] 2 \z J \ Z/ 

Here the O(z) terms vanish trivially in accord with our previous considerations, but we see from the Fierz 
identity that the 0(1) term is canceled as well. At four legs these cancellations continue to take place, but 
start to involve O(l) diagrams with a single quartic vertex insertion, for instance: 

^ 2 k ^2 k A2 

Evidently, such an intricate cancellation cannot be deduced by examining Feynman diagrams sep- 
arately. Let us look again at the 3-gluon vertex (3.5) from another point of view. Attaching a gluon 




propagator to the off-shell line obviously does not change the power of z, and the resulting off-shell 3-gluon 
current is a Lorentz vector. If we contract it with a simple fermion line, we obtain the first diagram in Fig. 
^, which is a part of a scattering amplitude — a gauge invariant object that is well established to behave 
as O (|) for the (k~,l~) shift. 
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Figure 3: Diagrams for the amplitude of 2 gluon and 1 fermion line. 

Moreover, the second diagram in Fig. |3| has zq momentum flow through its fermion line, so according 
to our previous discussion it is of order O (-) by itself. Thus we can conclude that the first one is O (-) as 
well. We obtained it by contracting the initial off-shell current vector with a correctly defined fermion line. 
The freedom of choosing the on-shell fermion momenta and helicities spans the whole Minkowski space. 
Therefore, the vector must have the same boundary behavior. 



k ■■■ I k ••• 1 




Figure 4: Contracting a gluon off-shell current with fermion lines. 

Along the same lines, we can now prove a very general statement: An n-gluon off-shell current with 
only two shifted like-helicity legs on shell behaves as O (|). The current has a free Lorentz index for each 
off-shell leg, so it is actually a tensor of rank (n — 2). If we contract every index with its own fermion line 
(independent of z), we will obtain an expression corresponding to a scalar amplitude (Fig. ||) which we 
know behaves as a whole as O Q) under the (k~ ,l~) shift. The freedom of choice of fermion momenta and 
helicities guarantees that if the contracted expression vanishes at z — > oo, then the initial tensor is bound 
to vanish too. 




Figure 5: Diagrams for 4 x {qq} g\§i with zq momentum flow through fermion propagators. 

Of course, we will still lack some diagrams to build the full amplitude, but the lacking terms will in fact 
be those for which we have already proven the good behavior at z — > oo. Indeed, the result resembles an 
amplitude for a process where (n — 2) (distinct) quark-antiquark pairs go to 2 gluons, so it should contain 
not only the diagrams which are given by the right-hand side of Fig. ^j, but also those where some fermion 
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lines have multiple fermion vertices and thus have fermion propagator insertions in them. Some of them 
look like the diagrams which are shown in Fig. || i.e. have zq flow through at least one of those fermion 
lines and thus vanish at z — > oo. 




Figure 6: Diagrams for 4 x {qq} — > gk<Ji, that can be reduced to the case of smaller number of gluon legs. 

Others, however, may look like the diagrams shown in Fig. ||, i. e. contain some fermion lines connected 
to the shifted gluons only through their connection to other fermion lines. These diagrams can be reduced 
to the case of a smaller number of off-shell legs in the initial gluon current. Thus, we can construct an 
inductive argument, for which we have already verified the base case of n = 3, to see that all the diagrams 
that we need to add to the contracted n-gluon current to form an amplitude behave as O (-) and the 
current itself is therefore bound to be O Q). 

By the way, this inductive proof did not use the weaker 0(1) statement of section 3.3.1, though we 
relied heavily on the O Q) statement of section 3.3.2. To conclude, let us recall the steps of our argument: 

1. Any diagram with zq momentum flow through at least one fermion propagator behaves well. 

2. The boundary behavior of the diagrams with zq momentum flow strictly through gluon propagators 
is the same as that of a gluon-only off-shell current. 

3. Any off-shell current with 3 gluon legs vanishes as z — > oo due to the cancellation which ensures the 
good behavior of the amplitude qq — > gg. 

4. Any off-shell current with n gluon legs vanishes as z — > oo to ensure the good behavior of the 
amplitude (n — 2) x {qq} — > gg, provided the good behavior of the (n — l)-gluon current and the 
diagrams with fermion propagator insertions. 

3.4 Mixed-helicity shift in special gauges 

For a generic gauge, the boundary behavior of individual Feynman diagrams under the [kl) shift is au- 
tomatically O Q) in the mixed-helicity case (k~,l£). But in the following section, we will find that in 
order to avoid unphysical poles we need to use special gauge choices = pi or n\ = Pk- Such gauges 
eliminate the z-dependence from one of the polarization vector denominators and thus turn the superficial 
behavior into 0(1). However, we can easily rephrase our power-counting arguments from the sections 3.3.1 
and 3.3.2 for the mixed-helicity case and find that the leading power of z must always involve at least 
one contraction of two of the following three vectors: zq x , [nklj^k) and (n^^lZ], with either g^ u or e\^ up . 
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Either by imposing n& = pi or n\ = pk, we can guarantee that any such contraction will give zero and thus 
ensure vanishing of the boundary term. 

It is worth noting that if we take = Pi and n/ = pk simultaneously, the superficial boundary behavior 
is worsened by two powers of z, and the argument will in general be invalid. Suppose that we first impose 
n i = Pk and have unfixed. Then we will be guaranteed to have no pole at infinity, but might still have 
an unphysical pole at Zk = [rikk}/\nkl}- Now if we take nk = pi, we can see that the pole Zk goes smoothly 
to infinity. In this way the unphysical pole and the boundary term can be traded one for another, and the 
problem is to find gauges in which neither survives. 

4. Avoiding unphysical poles 

In this section, we address the question of unphysical poles, i.e. the poles that come from polarization vectors 
( |3.1| ) instead of propagators. We construct explicit recursive proofs of the vanishing of the unphysical poles 
for the following currents: 

with n 3 = n 4 = • • • = n n ; 
with ?i4 = n-5 = • • • = n n = P3 ; 
With 713 = 71-5 = • • • = n n = P4 . 

It is straightforward to prove analogous statements for the currents with the opposite quark off shell 
or for flipped helicity assignments. For example, iJ(lg, 2q, 3~, . . . , (n— 1) , n+) has no unphysical poles 
under the [n— l|n)-shift if = ■■■ = n n _i = p n . One can also make a simultaneous flip of all gluon 
helicities trivially. 

In short, the good gauge choices are: 

• in the all-minus case, put all reference momenta equal to each other: m = q ; 

• in the one-plus cases, put reference momenta of negative- helicity gluons equal to the momentum of 
the positive-helicity gluon: n_ = p+. 

Note that in the one-plus case with the positive-helicity gluon in central position the unphysical poles 
vanish for a matrix- valued current \n~ . . . 5~4 + 3~| with both fermions off shell, i. e. lacking spinors on both 
sides of the quark line. In fact, there is strong evidence (see the numerical results in the following section) 
that it will continue to be true for such one-plus currents |n~ . . . (m+l)~ m + (m— 1)~ . . . 3~| irrespective of 
the position of the positive-helicity gluon with respect to the fermions. 

In each case, our recursive argument is based on the following expansion [1]. Consider constructing the 
n-particle current with one fermion line by attaching the n-th gluon to the corresponding (n— l)-particle 
current. Due to color-ordering, it can be coupled directly to the off-shell quark, to the (n— l)-th gluon, or 
to some more complicated gluon tree. If we focus our attention on those gluon trees that include the n-th 



1. [34)-shifted i J (l* Q , 2 Q , 3", 4", 5", . . . , nj) 

2. [43)-shifted iJ (l* Q , 2 Q , 3+ 4~, 5", . . . , nj) 

3. [34)-shifted ij(l^,2^,3-,4+,5-,...,n-) 
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gluon and attach to the quark line as a whole, we can expand the current according to the number of legs 
in such trees, as shown pictorially in equation (|4.l| ). 



4 5 ... 




4 ... 



71—1 



n-2 n-1 




(4.1) 



+ 



4.1 All-minus currents 

Let us prove that for the [34)-shifted all-minus current iJ (i*q,2q,3~ ,4~ ,5~ , . . . ,n~^j the residue at the 
unphysical pole z% = [n.33]/ [7134] vanishes when we make all the gluon reference momenta equal: 113 = = 
■■■ =n n = q. 

Now if we already know that the residue at the unphysical pole Z3 = [g3]/ [^4] vanishes for all the 
corresponding off-shell currents with fewer legs, then only the last two diagrams in fl4.1|) remain to be 
calculated. That can easily be done just by using color-ordered Feynman rules, where we make use of the 
Berends-Giele formula [1] for currents of like-helicity gluons, which in our conventions is given by 



iJ M (l~,2-,...,n- 



(4.2) 



\/%l][12]... [n-1 n][nq] 

Evaluating the sum at the pole 23 = 0, defined by [qo\ = 0, and performing some manipulations using a 
Schouten identity between the two contributions, we find the residue of (^4]) at the unphysical pole: 
4" 5- 4- 5- ... 




71 



+ 




i\q](m,n\q][<l\^2 ~ m)\2) 
(3|2|g][34][45] ... [n-1 n\[nq\ 



(4.3) 
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It becomes obvious that ( |4.3| ) vanishes due to the presence of the on-shell spinor |2) next to (2/2 — Tn). 

To conclude the proof, we do not even need to calculate the base of the recursion separately, because 
all the preceding formulas were general enough be valid for i J ( 1^, 2q, 3~, 4~ J as well. Indeed, in that 



case the last two diagrams in (4.1) turn out to be the usual Feynman diagrams with the Berends-Giele 
current representing just the polarization vector of the shifted 4th gluon. 



4.2 Currents with a single positive-helicity gluon in extreme position 



For the [43)-shifted current i J [Vq, 2q, 3j , A g , 5 g , . . . , n g j we put 



n r , 



P3 and rename n 3 = g. 

Consider the same expansion ( [4.1[ ). As in the previous case, for a recursive proof of the vanishing of the 
residue at Z3 = — (q3) / (q4) we only need to calculate the last two diagrams in ( fd.ip . We use the following 
formula for the one-plus Berends-Giele current: 



iJ^(l + ,2-,...,n- 



[l|7Ti,n|l] 



\/2[12][23] ... [n-l n][nl] 



^ PlPl^ + <21><lg>/' [4A) 



. 1=3 



in which we retain dependence on the reference momentum m = q of the positive-helicity gluon. This 
generalization is discussed in Appendix [B|. It turns out that relaxing one reference momentum results in 
only one extra term in (4^), which subsequently generates the pole for i/(3 , 4 , . . . ,n~) at (3q) = 0. 

Using the currents ([4. 2D and (|4.4D , we see again that the residue at 213 vanishes due to the presence of 
the on-shell spinor |2) next to ($2 — Tn): 




ii 



- 3 
+ 




-z|3](4 g )( g |f 3 ,n|3][3|(^ 2 -m)|2) 
(g|2|3](43)[34][45]...[n-l n][n3]' 



(4.5) 



This evaluation is valid as well for i J ( 1^, 2q, 3j, A g 



, thus ensuring the base of the recursive argument. 



4.3 Currents with a single positive-helicity gluon in next-to-extreme position 



In this section we consider a matrix-valued [34)-shifted current iJ ^1^, 2g, 3~, 4j, 5~, . . . , n~ \ with the 
positive-helicity gluon separated from the fermion line by one negative-helicity gluon. We put 77,3 = n.5 = 
■ • • = n n = pi and 724 = q. In expansion ( |4.1| ), we now need to examine the last three terms. More 
specifically, we need to examine only their g-dependent parts, since these are the ones that can have the 
unphysical pole at (g4) = 0. 

The very last diagram in ( |4.1[ ) contains the gluon current iJ (3~, 4 + , 5~, . . . , n~) , for which we do not 
know a simple analytic formula. Fortunately, according to the inductive argument outlined in Appendix 
[B], it does not depend on q, so that diagram cannot contribute to the residue at (g4) = 0. We are thus left 
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with the other two diagrams. As before, 
4 + 5" 



1* 



n 



\/2[43] 



i/(4 + , 5-,. 



.» ) / 2 "A +m J |4](3| + |3)[4| 



(P2 -P3) 2 - 



(4.6) 



We use the formula 



+ 



[n 3 3](n 4 4) ^ 

-m~n) + 



{p2-P3) 2 -m 2 



«3]<3| + |3)[n 3 | 



|4](n 4 | + |n 4 )[4| 

^(|4](n 4 | + |n 4 )[4|) - ^(|n 3 ](3| + |3)[n 3 



(4.7) 



2 [34] (43) 

which is derived simply from Feynman rules, to evaluate the third-to-last diagram in expansion (fOJ): 




n 



>/2(«4)[43] 



z\7(5 ,6 ,...,n ) 



(P2 - P3 - P4) 2 - m 2 



(4.8) 



|4]<?| + |9>[4| 



#2 ~ + m (q3) 
(P2 - P3) 2 - m 2 (43) 



|4](3| + |3)[4| 



After using the formulas (4.2) and (4^4) for the Berends-Giele currents and substituting 



la) = |4)M 

in the residue of the unphysical pole, we can combine the g-dependent terms of 
term with the following spinor matrix in the middle: 



(4.9) 

and (O) into one 



ii - % - #4 + m 



?4 + m 



+ 



#2 - % + m 



{P2 - % - Pa) 2 - m 2 {p2-%) 2 -m 2 (pa — P3 ~ Pi) 2 ~ m 2 (P2 - P3) 2 ~ m 



0. 



(4.10) 



Having established the induction, we turn back to iJ (l^, 2q, 3~, 4+, 5~J and see that its expansion 
(4.1) contains precisely the three diagrams that we have just examined in a more general case. This 



provides the base of our inductive argument for iJ [^-q, 2q, 3~, 4+, 5~, . . . , 

Our numerical results (see below) indicate that the statement about the matrix-valued one-plus cur- 
rents might be true irrespective of the position of the positive-helicity gluon as long as it is separated 
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from the fermion by at least one negative-helicity gluon. Unfortunately, it remains a challenge to show 
it. Here, we used an explicit formula for iJ ^1^, 2q, 3~, 4+^ to evaluate one of the diagrams. To prove 

the vanishing of the unphysical pole for iJ ^1^,2q,3~, . . . , (m— l)~,m+ , (m+l)~, . . . ,n~^j in the same 

manner, one would need to have an explicit formula either for iJ \1q,2q,3J , . . . , (m—l)g ,rnj \ or for 

iJ ( 1 0' 2 Q' m 9'( m+1 )fl '•••> r v)- 
5. Results for currents 

In this section, we apply the constructions established in the previous section to compute massive fermion 
currents from recursion relations. First, we list 3- and 4-point currents as a starting point. Next, we give 
a closed-form expression for currents with an arbitrary number of gluons if their helicities are all alike. A 
fully non-recursive version and its derivation are given in Appendix |C]. Finally, we state our numerical 
results for shifts producing recursion relations in the case of one gluon of opposite helicity. 

5.1 3-point and 4-point currents 

For completeness, we begin by listing the 3- and 4-point currents, which are straightforward to derive from 
Feynman rules, with full freedom of the choice of reference spinors. 



^(^^3;)=.™, (5.1a) 



iJ (%2 Q ,3-,4 ff -) = ? — i-^ * |4)K|2|3)[n 3 | - |n 4 ](4|l|n 3 ](3| 



[n 3 3][n 4 4] [ (3|2|3] 

+ m|4)[n 4 n 3 ]<3|+m|n 4 ]<43)[n 3 | j - | [n 4 3] ( |n 3 ](3| + |3)[n 3 | ) (5.2) 



+ [M(l'M](l| + |J}>h|)+^^(j/3-34)H|2)- 



Ull%,2* Q ,3-,4:- 



[34] 



[n 3 3][n 4 4] 
[n 3 n 4 ] 



(|n 4 ](4| + |4)[n 4 | 



n 3 ](3| + |3)[n 3 



(p2-P3) 2 -m 2 

ilk -fa) + [n 3 4] ( | n 4 ] (4| + 1 4) [n 4 1 ) + [n 4 3] ( | n 3 ] (3| + 1 3) [n 3 1 



(5.3) 
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Figure 7: BCFW derivation of i J 2 q , 3~, 4", . . . , 5 



[n 3 3](n 4 4> 



|l ](/; , K | <?j) ;i|) ^ ^ 2 +m 2 (|>^J(:i| + |:»[/»,| 
{P2-P3) ~rn z 



+ M(|4]<n 4 | + K>[4f) - ^) ( N(3| + |3}N 



2 [34] (43) 



[34] 



(43) 



(n 3 3)[n 4 4] 



[P2—P3)—m 



(5.4) 



(5.5) 



2(34) [43] vro ' ' (34) 
5.2 Closed form for all-minus currents 

For n > 5, we compute the all-minus current iJ (^q, 2q, 3~, 4~, . . . , by doing a [34) shift and setting all 
reference momenta equal to an arbitrary null vector q. Since we have established the absence of boundary 
terms and unphysical poles in the preceding sections, the BCFW expansion is given as follows: 



iJ UB' 2 Q' 3 ? ' 4 ? ' • • • ' n 9 = iJ[ l£, (2 - 3)q, 4 , 5 . . , n 



(h,h) 



n + m) 
9 I (p 2 -P3) 2 - m 2 

/J(l () .2 n ,3 3 ,(4 + 5)^6-,...,n- 



(P4 +P5) [ 



iM(-(2-3) Q ,2 Q ,3-) (5.6) 



; iM (-(4+5)^,4;, 5; 



fc=6 



See Fig. [5]. Because we are working with off-shell currents, the sum over intermediate gluon polarization 
states (h, h) must now include the unphysical polarization state combinations (L, T), (T, L), which vanished 
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automatically in the on-shell case due to the Ward identity. This subtlety was first treated in a similar 
context in [3]. Specifically, the numerator of the Feynman propagator is decomposed as 



(5.7) 



where 



k F 1 



It. 

kn 



(5.« 



It is clear that the gluon amplitude in the third line of (|0]) (the second line of Fig. 0) vanishes identi- 
cally, due to the form of the all-minus Berends-Giele current ( f4.2| ) contracted with any of the polarization 
vectors. The three-point gluon amplitude in the second line of ( j5.6|) (the last diagram of the first line of 
Fig. 0) vanishes as well. Therefore the only contribution that is left is the single term in the first line, 
involving a fermionic propagator. The general expression for an n-point all-minus current can then be 
written as 

\n~(n-l)~ . . . 4 _ 3~|2) =— — : ^ 

[q3}[qA]...[qn-2][qn-l)[qn] 



+ 



[n — 1 n] 



}(n\ + \n)[q\ 
[qn] ( |g](n|+|n 



n—l\ + \n—l)[q\ 



+ [gn-l](|g](n-l| + |n-l)[g| 



(5.9) 



{P2-Pz,n~2) 2 -rn 2 

^2-fzA + m 
(P2-P3A) 2 - m2 

rf2-$3+m 



(ra-2| + |ra-2)[g|) x ••• 
k/](4| + |4)[gf 



(p2-P3) 2 -m 2 
where the shifted momenta are defined recursively by 



M(3| + |3)M |2 



| A:] 



z k 

{kW 2 -F3,k-i\k + l} 
\k} = \k}-z k \k+l] 
\k+l) = \k+l)+z k \k) 

h,k = ?3,k - Z k (\k + l]{k\ + \k)[k+l\ ) , 



(5.10) 



with k = 3, 4, . . . , n — 2, and the initial values 

Z2 = 0, 



Z3 



(am 

(3|2|4] 



(5.11) 



- 18 - 



We have verified this formula numerically through n = 6 by comparison with sums of Feynman 
diagrams. 

The massless version (m = 0) was found in [1] for one helicity choice of the on-shell spinor, namely 
i J (^q, 2q, 3g , . . . , n~^j in our reversed fermion momentum convention, so that |2) = |2). In our calculation, 
rather than take the massless limit of ( |5.9|) , it would be more effective to return to the recursion relation 
as given in the nonvanishing first line of (|5.6| ), so that the propagators can be replaced by simple spinor 
products at each step of the recursion. Recovering the compact form of [1] is not immediate for general n, 
however, because we preserve a form of the current in which the quark spinor |2) is an explicit factor at 
the right of the expression, free to take either helicity value. This is important, because the shift of the 
quark momentum means that the full internal helicity sum occurs at each stage of our recursion. 

It is possible to solve the recursion exactly and write the shifted spinors for ( |5.9| ) in a fully closed form. 
We write this non-recursive form and outline its derivation in Appendix |C]. 

5.3 Numerical results 

Beyond the case of all gluons having the same helicity, we have found valid shifts numerically through n = 6 
in the case of one gluon of opposite helicity to the others (the "one-plus" case or its parity conjugate). 
A sufficient condition for a valid shift is to take the reference momenta of all the negative-helicity gluons 
equal to the momentum of the positive-helicity gluon: n_ = p + . 

For the choice of shifted gluons, we have identified two valid possibilities: 

• Shift the two gluons closest to the on-shell fermion; if they both have negative helicities, choose the 
shift so that the unphysical pole would come from the gluon adjacent to the on-shell fermion. (These 
shifts are all valid in the all-minus case as well.) 

• In the case with the plus-helicity gluon in central position shift the plus-gluon along with the any of 
the adjacent minus-gluons irrespective of their position with respect to the fermions. The unphysical 
poles then vanish, even with both fermions off-shell. 

To be more precise, we found that for iJ (]-q, 2q, 3~, 4+, 5~ , in gauge 723 = 715 = uq = P4 not 
only [34) shift produces no unphysical poles (as we have proved in Section 4.3), but [54) as well. Similarly, 
iJ (^*q, 2q, 3~, 4~, 5+, Qg ^ in gauge rc.3 = 724 = iiq = p$ suffers from no unphysical poles both under [65) 
and [45) shifts. 

In the 6-point case we also have currents with two plus and two minus helicities, but unfortunately we 
were unable to find a good gauge choice for them. 

6. Summary and discussion 

We have studied currents of n — 2 gluons of "mostly-minus" helicity and a massive quark-antiquark pair, 
where the antiquark is off shell. Because of the off-shellness of the antiquark, the choice of reference spinors 
plays an important role. 
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BCFW-type recursion relations are obtained under the following conditions, which ensure the absence 
of a boundary term and unphysical poles. The reference spinors of the negative- helicity gluons are all 
chosen to be equal. If there is a single positive-helicity gluon, its momentum is taken to be the reference 
spinor of the negative-helicity gluons. 

• In the case where all gluons have negative helicity, we have obtained both a recursive and a closed 
form for the current derived from recursion relations. 

• In the case where one gluon has positive helicity, and it is color- adjacent to the quark or antiquark, 
we have proven the validity of the recursion relation, but we do not have a closed form. 

• In the case where one gluon has positive helicity, and it is color- adjacent to two other gluons, we 
have found numerical evidence for the validity of the recursion relation in general, but were able to 
prove it only for the simplest configuration, with the positive gluon in next-to-extreme position. 

In Yang-Mills theory, an on-shell alternative to the BCFW construction is the MHV diagram expansion 
[27], in which maximally helicity violating (MHV) amplitudes play the role of interaction vertices, with 
a suitable on-shell prescription for the intermediate legs. For off-shell currents, there is apparently no 
sensible expansion in MHV diagrams when the off-shell leg carries color charge, such as the Berends-Giele 
currents for gluons. 

One might consider applying a BCFW shift to the massive fermion pair, but this construction fails 
off-shell. With a conventional definition of massive spinors [23, 26] in terms of a single reference vector, 
good boundary behavior is evident, but there are unavoidable, complicated unphysical poles, due to \fz- 
dependence of the denominators of the massive spinors. Even with both fermions on shell, the only shift 
known to be valid is quite specialized: each of the two massive fermion spinors has its reference vector 
constructed in terms of the other [19] . This choice is not well suited for repeated application in an analytic 
recursion relation, because it is undesirable to keep track of the data of internal legs. One would like the 
choice of reference spinor to be fixed once for all. Nevertheless, we looked at extending this construction 
off-shell. There is no z-dependence in the denominators, but when either of the on-shell massive spinors are 
stripped off, the miraculous cancellation reducing the boundary behavior from O(l) to 0(1/ z) no longer 
takes place. 

In the course of studying boundary behavior in Section 3, we have proven the good boundary behavior 
of general off-shell objects in Feynman gauge, as long as they contain at least two on-shell gluons that can 
be shifted. This meshes with a similar argument of [28] in the light-cone gauge q-A = specified by the 
BCFW-shift vector q, ( |2.3[ ). Thus we could see that it is not the boundary behaviour that hinders the 
BCFW recursion off shell, but the unphysical poles, coming from the polarization vectors. 

Several questions arise for future exploration. Is there any choice of shift and reference spinors that 
eliminates boundary terms and unphysical poles for more general helicity configurations? If so, can the 
recursion relation be solved neatly? Do some shifts give more compact results than others? Is there a 
neat solution for the current with a single gluon of opposite helicity, for which we have already proved the 



- 20 - 



existence of recursion relations? In cases where unphysical poles are present: is there any way to understand 
them, so that their residues could be incorporated explicitly in the recursion relation? Regarding the 
generalized Berends-Giele currents of Appendix |B|: can the current with one opposite-helicity gluon in a 
central position be written in a compact form, manifestly independent of the generic reference spinor? 
Further results addressing these questions would certainly illuminate our understanding of the BCFW 
construction and its applicability to gauge-dependent objects. 
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A. Ward identity argument in Feynman gauge 



In section 3.2, we used the Ward identity (3.3) to reduce the maximal superficial power of z at infinity by 
two, but we omitted some extra terms that are present in Feynman gauge. In this appendix, we fill in this 
gap in our argument. 

The supplementary terms we need to consider are due to the fact that the Noether current of the 
global gauge transformation receives additional contributions from the gauge-fixing and ghost parts of the 
effective Lagrangian. (Note that in the generalized axial gauge the gauge-fixing term tx (n^Afy 2 contains 
no derivatives and thus does not contribute to the Noether current.) At tree level, however, the ghost part 
of the current does not produce any non-vanishing diagrams, unless we consider a Green's function with 
external ghost legs, which is not the case. Thus the only non-trivial ingredient that we should worry about 
in Feynman gauge is the gauge- fixing contribution to the Noether current. 

To derive the Ward identities, we consider an infinitesimal global gauge transformation, 

A^A--gf abc a b A^ (A.l) 

which leaves invariant the -Rg-gauge-fixing Lagrangian: 

L ? = --L(VL^) 2 . (A.2) 

The latter generates the following contribution to the Noether gauge current: 

dL 



^ = g(Q^9f abc A c u = 9 -r hc A h »d v A™. (A.3) 



Taking a derivative gives 

<).. J"''i r\ - f 



d^{x) = 9 -f abc A b ^(x)d^d u A cu (x) 1 
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so we retrieve the following momentum-space operator: 

W t bcd f dA P 
(2tt)4 



Ut(i) = -ff b 



PvPp A cu {p)A d P{i-p). 



(A.4) 



This operator is to be inserted instead of the i-th leg and combined with the remaining (n — 1) legs of the 
off-shell current: 

k,fj,,a ... P,v,c 

l-p, p, d 



l _9_ rbcd 

' V J (27T) 



d4 P ^aur, 



I £ kp P" 




(A.5) 



Note that in ( |A.5| ) only the k-th leg is considered propagator-amputated, and we spell out its con- 
traction with the polarization vector explicitly. In the following we specialize to the Feynman gauge 
£ = 1 and once again neglect all color information. Since ( |A.4 ) has two gluon legs and already contains one 
power of g, then in order to construct a tree level contribution of order 0(g n ~ 2 ) from an object with {n + 1) 
external legs we need to contract two of them together without any interaction insertions. Thus the other 
(n — 1) legs must form a normal connected tree level diagram of order 0(g n ~ 3 ). An extra disconnected 
piece will naturally produce a J-function which will annihilate the integration in (A.5). 

We cannot contract together the two legs coming from (]A.4|) , because that would produce 5^(p + 1 — 
p) = 0. Moreover, if both directly contracted legs are not in ( [A.4] ), say the i-th and j-th, then we will 
get 8^\pi + pj) = 0. Obviously, we cannot contract directly a fermion with a gluon either. Thus we are 
left only with the options of connecting one of the two legs ( [A.4j ) with any of the remaining gluon legs - 
either the shifted one k or any of the unshifted legs. 

Contraction of the first leg of ( [A.4] ) with the on-shell gluon k vanishes immediately due to the transver- 
sality of the polarization vector whereas connecting the other leg to k results in p = k + / : 




(A.6) 



(2vr) V 4 ) (i-p + k)t k {k + i)„ ■ (k + iy 




k + l, v 



O 
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so that the only remaining element dependent on z is = O (j). Now if we take an arbitrary unshifted 
gluon leg with momentum pj and Lorentz index A and contract it with the leg p from (|A.4| ) , we get 




{2^\p+p 3 )-^p^ty J 



i+pj,p 




(A.7) 



In the second line we have written out the propagator of the (l — p) leg, making both z-dependent legs 
propagator-amputated, so that the diagram in the brackets is a standard O(z) and the overall expression 
obviously vanishes at infinity. 

The last remaining type of contribution is not so straightforward. Connecting an off-shell gluon leg to 
the (l — p) leg makes p equal to I + Pj and thus produces two more powers of z in the numerator: 




(A.8) 



(2vr)V 4 ) (f - p + Pj ) -^e (i + Pj y ■ t k {i + Pj y 



p 




l+Pj, v 



(l+Pj) 



Fortunately, what we can see on the right-hand side apart from other O(l) factors is just what we started 
with — a BCFW-shifted off-shell current with both shifted legs propagator-amputated and contracted with 
the polarization vector on one side and with the momentum (I +Pj) on the other. The main difference 
is that it has one gluon leg less than before and the momentum of the missing gluon is now added to the 
l-th. leg. 
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So we started with a superficially O(z) object with (n — 4) unshifted gluons on the left-hand side of 
(|3,3|), and we have managed to see that the right-hand side of the Ward identity gives us 0{\) contributions 
plus (n — 4) superficially O(z) objects of the same type, but with (n — 5) unshifted gluons. In the same 
fashion we can apply the Ward identity repeatedly until there are no unshifted gluons left, which proves 
that the initial object was indeed O(^). 

In this way we have verified that the Ward identity argument is still valid in Feynman gauge. 

B. Berends-Giele one-minus currents with a general reference spinor 

Berends and Giele [1] found a well-known recursion relation for multi-gluon currents. The explicit solutions 
they gave for currents with a single negative-helicity gluon require that it is color- adjacent to the off-shell 
leg, with a specific choice of reference momenta. In this appendix, we relax these conditions, for use in the 
proofs of Section 4. 

The Berends-Giele solution for the current i J^(l~ , 2 + , 3 + , . . . , n + ) assumed that the reference vectors 
rt2,...,n n for the positive- helicity gluons were all equal to the momentum pi, and that the reference 
vector n\ was equal to p2- Under these conditions, all the polarization vectors are proportional to Ai and 
hence mutually orthogonal, and the Berends-Giele recursion simplifies considerably. The orthogonality of 
polarization vectors is maintained even if we relax the second condition. Keeping n,2 = n% = ■ ■ ■ = n n = pi, 
but taking n± = q, we find the following compact result: 

v/2(12><23).->1) \{2l}{lq] £^ Pf.A-t J 

If the negative-helicity gluon is not color- adjacent to the off-shell line, we have a current of the form 
i J M (1 + , 2 + , . . . , m~ , . . . , n + ). If we choose reference vectors n\ = ni = ■ ■ ■ = n m -i = n m +i = ■ ■ ■ = n n = 
p m and n m = q, we find by an inductive argument that the current is independent of q. The base case with 
n = 3 can be evaluated explicitly as 

2- 3+) = ( 2 l^fi,3|2) [13] 2 
1 ' ' ' y/2{12) (23) [12] [23]' 1 ' ; 

As in the all-plus case, all polarization vectors are proportional to A m and hence orthogonal, giving the 
same simplified recursion (equation (5.47) of [1]), which does not contain the four-point vertex. According 
to the inductive assumption, it is enough to consider only the two contributions shown in Fig. || Using 
the result flETTD , we see that the ^-dependent terms from the two diagrams combine into a (/-independent 
contribution thanks to the Schouten identity. It would be interesting to find a compact form for this 
current, preferably one that is manifestly independent of q. 

C. All-minus formula 

In this appendix, we present the shifted spinors for the current with all negative-helicity gluons, equation 
( |5.9| ), in a fully closed, non-recursive form. The shifted momenta are written as follows, in terms of shift 
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1+ 2 + 




+ A 4 




Figure 8: Berends-Giele derivation of g-independence of gluon current with a positive-helicity gluon in central 
position. 

parameters Zk and new massless spinors \tk], both defined below. 

[k + l\k] 



The shift parameter is given by 



lk] = ltk] \k + i\t k y 

fe+l / k 

iw)=eij) 

i=3 \i=j 

f 3 ,fe = f 3 ,fc - « fc (|fc+l](fc| + |fc)[A;+l|). 



[t k \k + iy 



The spinor |^] is given by 

(k \ k / k \ 

Y[ Vj \(3\2\ + J2 n y*) K'l- 

j=4 / j=4 \i=j+l J 

Here we use the shorthand yk for the propagator denominator: 

yk = (P2 - -P3,fe) 2 - w, 2 ; 

the spinor [to^l takes the following values for k = 3,4, 

M=0, K| =m 2 (3|4|; 

for k > 5, 

where these elementary spinors [n^| and [vk\ are calculated as follows: 



[^1 



p=0 R 4 . fc - 2 
»l---»p 



(i P +i=k—l) 

p+1 



(C.l) 



(C.2) 
(C.3) 

(C.4) 

(C.5) 
(C.6) 



(-™y n %7(3i 2 i n ^-i+^-^+i) 

(<o=l) 



(C.7) 
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p=0 ^5,fc-2 



(i p+ l=fe-l) 
p+1 



-m 



2\P 



n n 

(io=3) 



(ij_l+2),ij ' P(ij+1) 



The sum is over all possible ordered pair lists ^ , defined by 



i? 



ii...i p 



where ii > i, i p < k, > ij + 2; and the complement sets 



(C.8) 



(C.9) 



(CIO) 



Finally, the products of /-matrices in ( |C.7| ) and must be ordered by the ascending numbering of the 

gluon momenta. 

C.l Derivation 

We now describe the derivation of the non-recursive formula ( |C . 1] ) , 

The expressions for \k] and \k + 1) follow directly from the recursions in ( |5.10| ), given the expression in 
flC.lD for Zk in terms of the new massless spinor \tk] • It remains to justify this expression and the definitions 
of \tk) , \u k ] , \v k ). 

Our derivation proceeds by induction on k. We take k = 5 as the base case. Here, even though the 
sets it! 4 ' 3 = S 4 ' 3 = 0, formula ( |C.7D already starts to produce a non-zero spinor [u^\ = — m 2 (3|2|^ z> 3 j 4|5| from 
the index p = 0, thanks to the boundary conditions io = 1 and i p+ \ = 4. We have verified with a direct 
computation that this value produces the correct expression in (|5.9| ). 

Assuming the validity of the formulas ( |C.2j ) - ( |C.§| ) for z^, [tk\ , [u^l , [vk\, we now demonstrate their 



validity for k — > k + 1. By the BCFW construction, in which only the first line of equation (5.6) survives, 
Zfc+i is obtained from z/. by making the following replacements of momenta, with the [34) shift: 

2 2-3 



1 3) 
1 3] 
4 
5 

k 



|4) 
|4] 
5 
6 

k + 1 



(C.11) 



We denote the results of the replacement ( C.ll ) on [tk\ ,[uk \ , [vk\ by [t' k \ , [u' k \ ,[v' k \. We now show that 
[ifc + i|, defined by 



fe +1 | = -(3|2|4][4| = -(3|2|4] 



'k+l \ k I fc+1 

n%- (4i2-3i+x n vAw\}' 

,j=5 I j=4 \i=j+2 



(C.12) 
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will once again satisfy the formulas ( |C.2| )-( |CL8| ), The prefactor — (3|2|4] conveniently reproduces the for- 
mulas for [uk\ and [vk\ and drops out of the ratio defining Zk+l- 

Applying the relation -(3|2|4](4|2-3| = y 4 (3|2| + m 2 (3|4| to the first term in QC.12D , we find 



•(3|2|4] [ J] yA <4|2-3| = JJ y 3 r ) <3|2| + JJ Vj ] [w 4 \. 

J=5 ) \j=4 / \j=5 



(C.13) 



The second term in ( p.l2| ) contains [w'j\ = [u'j\ + [v'j\. For the [v' k \ term, we have 



(3|2|4]K| = -m 4 E 



^j6,fc-l 



(i p+1 =k) 
p+1 



-m 



7 II Vi M-(3|2|4](4|5|} J] [JVi+D: 



»3 ' ^fe+1) 



(io=4) 



+™ 4 E E 



=0 R 6,fc-1 
^X-..»p 



(i p+ l=fc) 
p+1 



{-m 2 ) P J] W > <3|2|3+4|5| ]J [f fe _ 1+2 ),^( 



-m 



E E 

p=l ^j+fc-1 



J n-«P ) (io=4) 

)(i p+ l=fc) 
p+1 
<3|2| n [f fe _ 1+2)ji ,-y ( 



i=i 
(io=l) 



(C.14) 

where at the last step we included |3 + 4|5| into the ordered product, so the last sum goes over only the 
pair sets of the form Rf^^ 1 = (4, 5| U R ( ? ,k ~ 1 . 

A Cx...cp l ' j 1*2... Lp 



For the [u' k \ term, we start with 



"(3|2|4]KI 



[*?] 



-m 



E E 



(i p+ i=k) 
p+1 



-™ 2 Y II Vi {-(3|2|4](4|2-3|} [J [h-i- 



(C.15) 



fe+1) 



i=i 

(»o=2) 



where the hat over f > (i-_ 1 +2),i- indicates only that the ordered product of 's starts with |4 + 5 . . . |. We 

then apply the relation -(3|2|4](4|2-3|4+5+- ■ • | = y 4 (3|2|3+4+5 +• • • | + m 2 (3|4|5 H | to separate the 

terms in ( |C ■ 15| ) into two types: those starting with (3|2|3 + 4 + 5 . . . | and those starting with (3|4|5 . . . |. 
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Both can be easily incorporated into the ordered products just by shifting iq, giving 



■<3|2|4]K|=-m 2 £ 



p=0 p5,fe — 1 
"■il-ip 



[*?] 



(ip+i=k) 
p+1 



-m 



2\P 



n ^4 (3121 n [^-i+a),y-^+D 



■ _ 4,fc-l 
jG^... lp 



- 4 E E 



p=0 



3=1 
(<o=l) 

(i p+ i=fe) 
p+1 



(C.16) 



-m 



y n ^4 (3141 n 



•^(ij+1) 



._c,5,fe-l 



■4,fc-l m , , rr5,fc-l 



(io=3) 



where we have added {4} to ip = {4} U i by hand in order to include y± into the coefficient of 
the first term. The last term in ( |C.15[ ) is precisely 
Summing ( |C14 ) and ( C.16| ), we obtain 



-(3|2|4][u4| = [u k+ i\ + [vk+i\ = [w k+ i\- 
Finally, plugging (jCTTl and (|C17|) into (|Cl2|) we retrieve 

(fe+l \ fc+l / fe4-l \ 

n%- (3|2|+^ n w h, 

which is exactly our formula for I • 
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